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Abstract We consider a new approach to the problem 
of Bose-Einstein condensation (BEC) of polaritons for 
atom-field interaction under the strong coupling regime 
in the cavity. We investigate the dynamics of two macro- 
scopically populated polariton modes corresponding to 
the upper and lower branch energy states coupled via 
Kerr-like nonlinearity of atomic medium. We found out 
the dispersion relations for new type of collective exci- 
tations in the system under consideration. Various tem- 
poral regimes like linear (nonlinear) Josephson transi- 
tion and/or Rabi oscillations, macroscopic quantum self- 
trapping (MQST) dynamics for population imbalance 
of polariton modes are predicted. We also examine the 
switching properties for time-averaged population imbal- 
ance depending on initial conditions, effective nonlinear 
parameter of atomic medium and kinetic energy of low- 
branch polaritons. 

PACS 03.75.Lm; 71.36. +c, 42.50.Fx: 



1 Introduction 

Novel experiments on observation of Bose-Einstein con- 
densation (BEC) gives new opportunities to investigate 
of collective coherent phenomena in quantum and atomic 
optics and solid state physics - see e.g. [T], [3J. An im- 
portant tool for that can be taken with the help of the 
Josephson junction problem approach that has a long 
history relating to Cooper-pair tunneling in the super- 
conductor Josephson junction (SJJ) systems - see e.g. [3]. 
Various aspects of Josephson phenomenon and super- 
fluid behavior for two weakly linked atomic condensates 
in double-well potential have been discussed in a num- 
ber of papers @]-[7|. However just recently the coher- 
ent oscillations and self-trapping regime have been ob- 
served experimentally in [5] for two condensates of 87 Rb 
atoms placed in a special optical dipole trap. The prin- 
ciple item is that photonic "condensate" exhibits a simi- 
lar nonlinear behavior in the problem of optical Joseph- 



son junction (OJJ) consisting of two tunnel- coupled op- 
tical fibers and/or semiconductor waveguides (quantum 
wells) [5] . Another type of Josephson coupling (nonlinear 
Josephson effect) for condensate modes has been con- 
sidered in [10J. In particular, nonlinear dynamics like 
macroscopic self-trapping and intermodal population ex- 
change, has been theoretically analyzed for condensate 
collective modes (eigenmodes) , and namely, for ground 
state and first excited state modes. The effective cou- 
pling of these modes is guarantied by atom-atom scat- 
tering, i.e. by the own nonlinearity of BEC placed in 
highly anisotropic trap. 

However, the observation of mentioned above phe- 
nomena with pure atomic BEC medium is strictly lim- 
ited by extremely low temperature (about hundred nK) 
conditions. At present the problem of high-temperature 
quasi-condensation phenomenon that occurs for cavity 
excitons and polaritons in so-called strong coupling regi- 
me evokes a great interest p], [TT|-[T3]. Speaking more 
precisely the low branch cavity polaritons can be consid- 
ered as a weakly interacting two-dimensional Bose-gas 
for which the Kosterlitz-Thouless phase transition oc- 
curs at high enough temperatures due to the extremely 
small mass of polaritons m e g = 5 • 10~ 33 g under con- 
dition of strong coupling regime for excitons and quan- 
tized e.m. field. Although the polariton condensate is 
the non-equilibrium state and occurs in a specially con- 
trol situation, the macroscopic occupation of ground 
state with momentum k = has been examined experi- 
mentally for quantum wells inserted in a semiconductor 
(CdTe/CdMgTe) microcavity - see [T2J and more recent 
results in [13| . The coherence properties have been ob- 
tained in [11] for low-branch polaritons above the thresh- 
old of generation in semiconductor (Ga-As) microcavity. 
It has been demonstrated that statistical properties of 
the excitons can be converted to match the excitons in 
coherent states. 

The problem is that the BEC state in semiconduc- 
tor systems, fabricated at present, is essentially limited 
by disorder factor. Recently it was suggested to use a 
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strong disorder in the microcavity to obtain polariton 
quasi-condensation |15j . However it is not clear until now h( 
it is possible to achieve true condensation for exciton- 
polariton substance in semiconductor microcavity for this 
case. 

In our paper [16j it has been shown that polariton 
true condensation can be achieved for relatively hot atomic 
vapor in the case of polariton trapping, when atoms 
are strongly coupled with optical field in the resonator. 
The degeneration of weakly interacting 2D-gas of po- 
laritons can be also associated with the effect of so- 
called "spectrum condensation" observed previously for 
sodium atoms in the cavity. 

For some practical purposes the coherent properties 
of upper branch of polaritons can be important as well. 
In fact, the problem of formation of quantum entan- 
gled polariton states from different branches [T7] is the 
case in respect to construct the algorithms of quantum 
cloning and storage for optical field into the cavity po- 
laritons being recently proposed [18] . However, in the 
experiment [14] the upper polariton branch has not been 
developed in the photoluminescence when the quasi-con- 
densation of lower polariton branch takes place under the 
thermal equilibrium condition. The existence of large en- 
ergy gap (about ten of meV) between upper and lower 
branches is also a reason for that. 

But the macroscopic population of the ground (low 
branch) polaritonic state in atomic system can be only 
one of possible way to obtain interbranch dynamical 
regimes for such a polariton system. In particular, the 
energy gap in atomic system can be in 10 4 — 10 5 times 
smaller then for semiconductor one [16] . [19j . and some 
(small) population of upper branch does exist even un- 
der the quasi-equilibrium conditions. Upper polaritonic 
branch can be populated also by using Landau-Zener ef- 
fect when the tunnelling phenomenon between two eigen- 
states of the system takes place for the time-dependent 
detuning - e.g. [2T] , 

Another possibility to observe such a dynamics is de- 
termined by a special scheme of coherent excitation of 
polaritons from different branches [22] . Quantum beats 
and Rabi oscillations in exciton-photon system have been 
observed in this case |23j . 

In present paper we consider the problem of nonlin- 
ear dynamics for polaritons of different branches under 
the mean field approximation. Our approach has some 
analogy with the problem of macroscopic quantum tun- 
neling and interference that takes place for atomic con- 
densate being associated with ground state and first ex- 
cited state. Actually, low and upper branches of polari- 
tons represent two types of normal modes in atom-field 
interaction. They can be coupled to each other due to 
atomic nonlinearity. 

The paper is arranged as follows. In Section 2 we es- 
tablish nonlinear model of interaction of two-level atoms 
with quantized e.m. field in the cavity for both momen- 
tum and coordinate representation - cf. pQ, [19], [24] , 



The adiabatic approximation is discussed for the last 
case. In Section 3 we propose a new approach to de- 
scribe a nonlinear Josephson junction problem for both 
upper and lower macroscopic polariton modes coupled 
in momentum space due to nonlinear (Kerr-like) atom- 
field interaction. The analysis of Josephson dynamics 
and switching of various regimes of polariton interac- 
tion is presented in Section 4. In Conclusion we briefly 
summarize our results. 

2 The nonlinear model of atom-field interaction 
in the cavity 

To describe the interaction of two-level atomic system 
with a quantum e.m. field, we consider a localized exci- 
ton model within the framework of the Dicke-like Hamil- 
tonian for the interaction of the N two-level atoms with 
a quantized e.m. field (cf. pQ, [2], [19]): 

H = J2E ph (k) + (b]bj ~ otoj) 

k 3 




where ipk y^kj ^ s ^ ne annihilation (creation) operator for 

photon with momentum k; a (a^) and b (&^) are the an- 
nihilation (creation) operators for atoms at the lower \a) 
and upper |6) levels, respectively, N is the total number 
of particles, E at is the energy of atomic transition be- 
tween the | a) and \b) levels (here we neglect the motion 
of atoms in the cavity), -E p h (k) defines the dispersion 
relation for photons in the cavity, and g characterizes 
the atom-field coupling. The last term with coefficient 
k in the Hamiltonian JT]) takes into account a nonlin- 
ear atom-field interaction, and can be considered in the 
framework of nonlinear Jeans-Cummings model. The in- 
teraction implies a polariton-polariton scattering from 
different branches, and is determined by Kerr-like (opti- 
cal) nonlinearity of atomic system. 

The character of interbranch polariton interaction is 
determined by sign of effective parameter g ■ kN / . 
In particular, repulsive particle interaction (with positive 
scattering length) is realized when g ■ kN < (E at > 0). 
On the other hand, attractive interaction (with negative 
scattering length) occurs for g ■ kN > 0. In the paper 
we take into account the positive coupling parameter g 
only. 

In contrast with semiconductor systems we neglect in 
Hamiltonian JT]) by terms describing the nonlinear pro- 
cesses of exciton-exciton scattering - see e.g. [Tj. Such 
processes can be reduced to parametric amplification of 
polaritons and/or stimulated emission - see, e.g. [2~1] . 
Our approach is obviously valid for ensemble of non- 
(or weakly) interacting atomic gas for which the average 
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interatomic separation length is much larger than char- 
acteristic length of interaction between the atoms [25] . 

For high-reflection mirrors in the cavity the normal 
(to the plane of the mirrors) component of the photon 
wave vector k± is quantized, i.e. fcj_ = irm j X cav , where 
L cav is the cavity length, m is the number of modes. For 
the mode component in parallel to the mirror plane fcy we 
have a continuum. In paraxial approximation (fcy <C k±) 
the dispersion relation for photon energy i? p h (k) inside 
the cavity can be represented as follows: 



2 U2 



Eph (k) = he |k| = he 



ki 



k 2 w hck i 



h z k 



2m 



P ii 



(2) 



where m pn = hk±/ c is the effective photon mass in the 
cavity. 

In the paper we impose a strong coupling regime for 
which the cavity modes are coherently coupled with col- 
lective atomic exitations that can be introduced with the 
help of annihilation (<fi) and creation (cfy) operators as 

(cf. my- 
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j=i 
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b)a 3 



N 
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(3) 



The low (or zero) density limit which is used for ex- 
citons in semi-conductor structures (see e.g. PQ) corre- 
sponds to two-level atomic system without inversion. In 
this case the atomic excitation operators <j> and cfy in 
Eqs. © obey to usual commutation relation [</>, (jy \ ~ 1 
for Bose system when the atoms occupy mainly a ground 
state | a) of energy. 

Let us consider the system ([I])-® of atom-field inter- 
action in the cavity in coordinate representation. In this 
case we can introduce appropriate Lagrangian density L 
in the form: 



L = 
U 



U 
ih 



L at +L 
(yd** 

*2 



int j 

qr* 



~&t) 

oz oz ) 



2m pn 
L at = ^ ( 

L int = .9 (^*<£ 



(4) 



1*1 



where the & = & (r, t) function (i.e. the order parameter) 
describes both the time evolution and spatial distribu- 
tion of optical field in the cavity. In paraxial approxi- 
mation ((2]) under the mean field theory condition the 
value of ^ = <P (t) is the C-number that characterizes a 
macroscopic excitation in time of atomic system (here 
we neglect the motion and space distribution of atomic 
cloud in the cavity). In expression ^ we also introduce 
the external (harmonic) potential t/ pt = C^opt (x, y) for 
the photon trapping - cf . [TB] . 

Using the expression ^ under the variational ap- 
proach (see e.g. [9]) we obtain the following coupled 



equations for the field ^ and the atomic excitation $ 
respectively: 



ih 



d 
M 



'Wz 



'2.m p h 
-g$ + k (2 I 



<f = 

d 2 
dx 2 

\Sf\ 2 $. 



dy' 



+ U opt (x,y))V (5) 



ih— = E at <P 
at 



■g& + n\W\ V 



(6) 



Let us consider a solution of Eqs. ([5]), (HJ) under the 
adiabatic approximation for which we suppose d<P/ dt = 
for atomic excitation. Eliminating the <P variable from 
Eq. ([6]) we arrive to self-consistent generalized Schrodingcr 
equation for the order parameter W that has the form: 



ih 



d 



+U opt (x,y) 



u W\ 



h z 

2m p h 
'' + dW\ 



d 2 
dx 2 

9. 



d 2 
dy 2 



(7) 



where u — —Agn/E^, d — — 2>n 2 / E at . Here we also omit 
the term with g 2 /E at that can be also eliminated with 
help of simple unitary transformation. 

Equation ([7]) describes a behavior of optical beam 
interacting with nonlinear atomic medium in the cavity. 
The last term in Eq. can be discarded within the limit 
\d\ \Hf\ <§; |u| that corresponds to phenomenological 
(macroscopic) approach for weakly interacting Bose-gas 
of photons in the cavity - cf . . The coefficient u that 
is determined by Kerr-like nonlincarity of atomic vapor 
is a control parameter for that. 

The coefficient d that is always negative in Eq. ([7|) 
characterizes a quintic nonlinearity of the medium. Re- 
cently a superfluid behavior has been considered for such 
a system in [27]. In particular, formation of stable vortex 
lattices for "photonic condensate" has been numerically 
proposed even for the case without the external trapping 
potential U op t (x, y) use. 

In general case when both relaxation rate of optical 
field If and atomic excitations dephasing rate r cx are 
under consideration the adiabatic approach is valid when 
an inequality 



r ox » if 



(8) 



is hold - cf. [IS]. 

In the present paper we consider the atom-field in- 
teraction in the limit when the rates r cx and If are the 
same order of magnitude. In this case we can examine the 
system described by Hamiltonian (TT]) using in fact the 
two macroscopically populated polaritonic modes from 
different brunches. 
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3 Boson Josephson junction (BJJ) model for 
polaritons in the cavity 

Let us consider a Josephson dynamics for two modes, 
i.e. upper and low branch polaritons, those are coher- 
ently excited in the cavity - cf. [33J . We suppose that one 
quantum optical mode (with k = as a limit) is only 
macroscopically occupied. Thus, in Hamiltonian ([I]) we 
switch over to macroscopic (continuous) variable ip (k) 
for optical field describing this mode. In this case we 
introduce the time- dependent classical wave functions 
*i ja (t) = yJMt)e' 9 ^W 



as: 



1 



#! = -j= (V(k) - 0) 

<2> 2 = ±= (V(k) + 0) 



(9) 
(10) 



that characterizes the bosonic quasi-particles in atomic 
medium corresponding to upper and lower branch po- 
laritons; Ni t 2 (t) is the time dependent average number 
of polaritons for upper (iVi) and lower (N 2 ) branches 
respectively, 9\ and 9 2 are the phases. 

Now we introduce the new variables <r and 9 as 



N!-N 2 



(11) 



that characterize a fractional state population imbalance 
(s) and relative phase (0) for polaritonic modes respec- 
tively, N cx = Ni + N 2 is the total number of polaritons 
in the cavity. 

In the paper we pay our attention to the problem of 
polariton interaction along the cavity length under the 
atom-field resonance condition A = E ai — %c\k±\ = 0. 
The Hamiltonian ([I]) , describing the superfluid dynamics 
of the polaritons in momentum representation, can be 
evaluated in the terms of new variables <; and 9 as 



H 



: A cS fc 2 ~ A cS (3 (Vl - ^COS# - y/TTW) 

x\/l - c, 2 cos9 + AE cS c,, 



(12) 



where (3 = E tI / \g\ is the scaled (normalized) energy 
of polaritons in transverse plane, E tI = fi 2 ^ 2 ^ (2m po i) 
is their kinetic energy producing the effect of quantum 
pressure, m po i ~ 2m p h is the mass of polaritons in the 
cavity. It is possible to define Josephson-like length Aj = 
\f\P\ / fcy that specifies a spatial control parameter in the 
system - cf . [3J , [5] . For strong coupling regime with the 
value of atom-field coupling parameter g ~ 0.05 meV 
the length Aj ~ 11 /im. The result is the same as for 
atomic BEC in the BJJ problem - cf. [5]. The crucial 
dimensionless parameter A c g = nN ex / [2\g\ (l + /? 2 )] = 
A / (l + P 2 ) determines a nonlinear Josephson coupling 
between polaritons of different branches (the parameter 
A < 0, i.e. «iVe X < 0, for positive scattering length of 
polaritons). 



In expression (| 1 2|) for Hamiltonian H we define an 
effective zero-point energy difference AE e g as 



AE eS = a- ylofr v 7 ! - 7 2 (l - P 2 ) cos 9, 



(13) 



where we have made a denotation a = yl + (3 2 (1 + yl e ff)- 
Physically, the parameter a characterizes the energy gap 
between polariton branches, i.e. determines the unalter- 
able zero point energy difference when condition /3 = 1 
takes place. 

Using the Hamiltonian ([T]) and definition (fT"3"|) for 
AE c s one can obtain the following equations for canon- 
ically conjugated variables ? and 9: 



§ = -A cS (? (1 - p 2 ) - P^JT+W 

+2/V1 - ^ 2 cos 9) \/l-<r 2 sin 9, 



(14) 



d0 
7h 



A r 



+2A cS p ? 



1 - 2q 2 



cos 9 



Vl + f3 2 s 



scos 2 9 (15) 



l-2^ 2 



cos 9 



where r = gt/Ti is the dimensionless time coordinate. 

Thus, the set of equations (fT4I) . p3|) represents an 
extension of the BJJ model in nonlinear dynamics ap- 
proach with the phenomenon of tunneling for polari- 
tons from different branches. In the case P = the 
Eqs. (fT4l) . (fT5]) describe the polariton arising effect along 
the cavity length only. 

Obviously, it is possible to establish some link be- 
tween the system of Eqs. (TT41 . (fT5")) describing the po- 
lariton dynamics in the cavity and the other Josephson 
dynamical systems - cf. (3j - [10] . In particular, the equa- 
tion p4p characterizes the interbranch current / for po- 
laritons. In general it has a complex form. However we 
can represent the current I for small values of P as 



-/cvT 

+ L I 



1 



2 sin6» + Ji<rVT 
V) sin 2(9, 



C sin I 



(16) 



where I c — \n\N eK E% T / (2 \g\ K) is the positive critical 
current for polaritons in transverse plane (we took A e R < 
0), Ji = |k|AT cx / (2Ti) is the current characterizing the 
polaritons along the cavity length. 

In general the interbranch current I is essentially 
nonlinear one; the case is quite different from the Cooper- 
pair current in the SJJ problem - cf. [3J, [4]. From Eq. 
(I16|) it is easy to see that the first term, describing the 
supercurrent of polaritons in transverse plane, is com- 
pletely similar to the BJJ inertrap current occurring for 
atomic condensates. The last term in Eq. (| 16|) represents 
a nondissipative current determined by interaction of the 
polaritons inside the cavity. Physically, this term can be 
associated with the second order contribution to Joseph- 
son current due to the effect of tunneling for atoms in 
condensate state to noncondensate quasiparticles in the 
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BJJ problem [7]. For small population imbalance |<r| <C 1 
the contribution to the current of polaritons along the 
cavity length can be neglected . 

Energy difference AE c g defined in Eq. (flU)) plays a 
principal role in total polariton dynamics. In respect of 
the atomic BJJ problem the zero-point energy differ- 
ence AE c ff is determined by both the geometry of the 
trap and the nonlinearity of the atomic condensates [5] , 
[Hj. The energy AE c g, determined by applied voltage 
across the junctions, results in oscillations being in anal- 
ogy with the ac-Josephson frequency w ac ~ AE e g in 
the SJJ problem. However in the case of polaritons, the 
energy difference AE c r in Eq. (fT3| is mostly depends 
on dynamical variables q and 9 in general. Physically, 
such a behavior is caused by nonlinear (Kerr-like) cross- 
interaction of polaritons from different branches in the 
cavity. The energy difference AE C & defined in Eq. (fT3|) 
can be suppressed for the case A e g < 0. For example, we 
have AE e s = for the values p = 1, A c g = — 1 what 
implies the polaritons with momentum fcii ~ 1 /Aj and 
kinetic energy E tr = \k\N cx / A . This condition is com- 
patible with the requirement for definition of the heal- 
ing (or coherent) length £ c that plays an important role 
for superfluid systems - cf. [25| . In particular, a typical 
length scale £ c for that is determined by the condition 
when the E iT energy is of the order of the interaction 
energy for condensate atoms (determined by the reiV ex 
parameter). For the case under discussion we have the 
coherent length £ c ~ Aj. 

4 Nonlinear superfluid dynamics for trapped 
polaritons 

4-1 Small- amplitude oscillations 

We start our analysis from the case of non-interacting 
(ideal) gas of polaritons when A — 0. The solutions of 
Eqs.JHJ), (fT5j) are simply presented in the form: 



6 = g^+P^t = —t, C = ?({)), (17) 

where AE = 2gy / l + P 2 is the energy splitting for non- 
interacting polaritons, ^(0) is the initial population dif- 
ference. Thus, in the case under consideration the po- 
lariton branches are not overlapped and therefore the 
polaritons arising in the problem are still uncoupled. 

Now let us consider a small- amplitude oscillations 
case for population imbalance q when the time average 
values (9) = ("zero phase" oscillations) or (9) = n 
("7r-phase" oscillations). Linearizing the Eqs. (|14p . (fT5)) 
in respect of the q and 9 variables we arrive to equation 

+ Ql v q = iV , (18) 

where two parameters, i.e. dimensionless angular fre- 
quencies of oscillations i?o,7r and external forces Fq )W , 



are defined as 

% = (4 ,±4 ,) 1/2 , (19) 

r? JP o,. = y/ (20) 
o,* = VA* (l-P 2 )(aTA cS ±A cS p 2 ), (21) 
F 0;7r = -Ef^ (a T A cS ± A cS p 2 ) . (22) 

In expressions (fT9|) -(|22 |) we introduce two characteristic 
(dimensionless) parameters, i.e. the Josephson coupling 

energy #f ,7r) = A cfi p (2 T \A + P 2 ) - h, and "capaci- 

tive" energy £ , ^°' 7r) = A cS p (a =f \/l + P 2 ^j that describe 
the polariton nonlinear interaction for different branches 
- cf. i. 

The expressions (fT9]) -(l2~T j) determine the dispersion 
relations for a new type of collective excitations arising 
in the cavity for polariton interaction (we stress here the 
fact that /3-parameter is proportional to momentum as 
P ~ fcjj). Nevertheless it is useful to establish some link 
between the problem considered here and some other Jo- 
sephson junction systems. In particular, the small ampli- 
tude oscillation frequency regime, described by Eq. (fl9|) . 
has the same structure as a frequency of small oscilla- 
tions for the atomic population imbalance occurring for 
Bose-Einstein condensate confined in the W-potential :5j . 
The second term in Eq. I|19|) characterizes the linearized 
Rabi oscillations with frequency X?r, o,tt of the population 
imbalance. The angular frequency I7jp in Eqs. (fl9|) . ([20]) 
can be considered as an analog of the plasma frequency 
for Josephson plasmon in the problem of a Cooper pair 
macroscopic tunneling process across the junction - cf. [3] , 
|29j . The driving force Fq^ depends on /3-parameter that 
determines quantum pressure of polaritons in transverse 
spatial plane. 

We distinguish two principal regimes of collective am- 
plitude oscillations for the cavity polaritons. Namely, 
within the limit 3> i?jp so called Rabi regime oc- 
curs. For quasiparticles with a small momentum fcii (i.e. 
for p 4 <C 1) the population imbalance oscillates with 
angular frequency i? ps ^A.5Al s P 2 — 0.5A c gP 2 + A cS 
vs r under the zero phase regime. In the case of ab- 
sence of polaritons in transverse plane (/3 = 0) we have 
J?o = %o = VA (for A > 0). The expression ([20]) 
and numerical simulation show that the Josephson cou- 
pling energy Ej is vanished, and the contribution of Rabi 
frequency J?r in Eq. (|19p is important for parameter 
/3-> 1.7. 

On the other hand, the Josephson regime is realized 
for the polaritons in the cavity within the limit -C 
f?jp. Such a regime of polariton oscillations is achieved 
for the parameter p = 1. 

In the limit of large P (P 2 ^> 1), i.e. for the 
quasiparticles with large momentum fcn , the frequency 
of amplitude oscillations is Qq ~ y/ ' A 2 S P A + \A c s \ P 3 = 

\J \J\A cS \ PE C Q K + (for A cS < 0) that is in good 



6 



A. P. Alodjants et al. 



agreement with the oscillation frequency expression ob- 
tained for atomic population imbalance in the BJJ prob- 
lem [29]. Thus, crossover to Rabi oscillations with fre- 
quency X2r, ~ y/]Aes\W (f° r \A e g I "C 1) can be easily 
achieved by varying, for example, of the /1-parameter 
(i.e. the intensity of laser field). 

The solution of Eq. (I18|) is given by following expres- 
sion: 



S0,tt = C COS ( A),7r"r) + <J 



(0,7r) 

dis ' 



(23) 



where C is the constant determined by initial condi- 
tions, (<;o,7r (t)) = sjjfa = -Fo,tt /^o.tt characterizes the 
displacement for the time average fractional population 
imbalance (<;) in respect of the zero value. The displace- 
ment ^ig 7 ^ = and the driving forces Fj — for (3 — 
or = 1, A c g — —1. In this case small amplitude oscil- 
lations of the ^ variable occur around the value (<;) = 0. 

Developed above approach can be violated, for ex- 
ample, for running-phase regime when, first, the phase 
difference 9 grows unrestrictedly and, second, the essen- 
tially nonlinear oscillations of population imbalance ? 
take place. We are also emphasis here some difficulties 
to obtain so-called Fock regime by using the Hamilto- 
nian ([I]) when fulfillment of the inequality E c ^> Ej is 
required - cf . [29] . 



4-2 Stationary states of polaritons 

First, we focus our attention on the out-of-phase station- 
ary states of polaritons in the cavity determined by the 
conditions sin (9 ^ 0, cos# ^ 0. From (fl"4|) . (fl~5l) it is easy 
to obtain the following solutions 



A + 2 

A VW 2 



1 



arccos 



A - (3 2 + 1 



(24) 



sgn (A) y^ 2 - 4/3 2 [A + 1) ' 



that represents unstable saddle point with energy Hq = 
— ^ (A + 1). The existence of solutions (f2"4"]l implies ful- 
fillment of conditions for population imbalance <r and 
phase 8, i.e. \s\ < 1, |cos#| < 1. The fact results in cer- 
tain condition for parameters /3 and A that looks like: 



A<-\(f + l) 



(25) 



Second, we consider another out-of-phase solution of 
Eqs. (I14p . (fT5"]) with sm.9 — ±1. For stationary fractional 
population imbalance <^ and energy H of the system we 
find: 



a,2 = ±1, 
2 + l 



±2/3 + V/3 2 + l V y/P + l 



1± 







(26) 



_ /?V^TT _ {A + (3 2 + l) y/0 2 + l 

' (27) 



„ _ A 



The solutions (f2"o) . (j2"7) do exist if the parameters /3 
and /l fulfill the conditions : 



A 



(0 2 + l) 



3/2 



±2/3 + V/3 2 + 1 



for ft, 2, 



and 



= /3 2 - 1 for ftj 



(28) 



(29) 



where /3 € 



Thus, the solutions ([H, $27$ 



l . l 

'7S' vlj 

represent the unstable saddle points in phase portrait. 

Third, Eqs. ([M]) , (fl~5|) have zero-phase (in-phase) and 
n-phase stationary solutions with the phase 9 deter- 
mined by the condition sin (9 = 0. The stationary points 
in this case can be found out from solution of algebraic 
equation 



A,, 4 



A 2 e 



A 3 s 2 + A^ + A 5 = 



(30) 



for fractional population imbalance where the real co- 
efficients Aj = Aj (A, (3) are the functions of A and /3 
(we are not directly represent them in the paper due to 
complexity of the relevant expressions) . In general, there 
exist four roots of Eq. (|3"U|) of which only the real solu- 
tions are of interest. Note, that one of the solutions of 
(j3TJ]) may become a saddle point if the condition (|2"5j) is 
not fulfilled. 



4-3 Phase plane portrait description 

From mathematical point of view the equations (Q3 
(|15p describe an undamped nonrigid pendulum with phase 
9 and canonical momentum pg = <;. In the paper we 
consider a phase-space description to understand the dy- 
namics of the polariton interaction in the cavity. In Fig.l 
the evolution of the phase portraits for various values 
of normalized kinetic energy f3 of polaritons is presented 
for the same nonlinear yl-parameter. The trajectories are 
plotted for different values of initial conditions for the 
population imbalance and the phase with A and (3 pa- 
rameters kept constant. Closed loop trajectories in the 
gray area marked in Fig.l correspond to periodic oscil- 
lations of pendulum coordinate and phase in terms of 
nonlinear pendulum analogy. The stationary points Sj 
(j = 1, ...,6) correspond to the roots of Eq. (j30|) . White 
color regions 1 in Fig.l determine the running phase 
regimes that correspond to anharmonic rotation of non- 
rigid pendulum with relative phase —oo < 9 < +oo. Sta- 
tionary saddle points Pi. 2 at the separatrixes in Fig. la, b 
are determined from expressions (|24p respectively. 

For arbitrary value of zero-point energy difference 
AE c fi defined in (|13[) the phase-plane portraits are asym- 
metric as a rule. However, the phase portrait established 
in Fig. lb becomes completely symmetric for AE c g = 0. 
In particular, the small amplitude zero phase Joseph- 
son oscillations of fractional population imbalance take 
place here in the region 5 with the frequency fi^ = i?jp 
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Fig. 1 Phase-plane portraits of the variables ? and 6 /it for: 
a) (3 = 0.15, b) (3 = 1, c) (3 = 2.08475. For all cases the 
parameter /l = —2. The points Sj and Pj are the station- 
ary points. The regions numerated as '1,2,. . . ,5' determine 
the various dynamical regimes separated by dashed curves 
(separatrixes). 



around the average value (s) = with zero displace- 
ment. The fact is in agreement with Eqs. (fl"9|) - (|23|) . Two 
symmetric regions numerated as 2 in Fig. lb character- 
ize the large-amplitude (7r-phase) oscillations of pendu- 
lum around the top of vertical axis. A complete energy 
transfer between polaritons from different branches takes 
place in this case approximately. For small initial val- 
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T 



Fig. 2 Fractional population imbalance q against scaled 
dimensionless time r. The parameters are: /3 = 2.08475, 
?(0) = -0.97, 6»(0) = for curve 1, p = 0.15, c(0) = -0.5, 
0(0) = 7T for curve 2, /3 = 2.08475, ?(0) = 0.9, 0(0) = vr for 
curve 3 and (3 = 1, ?(0) = -0.259, 0(0) = for curve 4. For 
all cases the parameter A = — 2. 

ues of the fractional population imbalance <; 2 (0) <C 1 
we obtain the small-amplitude 7r-phase oscillations with 
frequency f2„ - see (fl"9|). 

Other regions in Fig. lb and all of them in Fig. la, c 
specify so-called macroscopic quantum self-trapping regi- 
mes (MQST) for which a fractional population imbal- 
ance <r oscillates around the non-zero value (<;) 7^ 0. It 
implies that population of either of polariton branches is 
never completely depleted. In fact, in our case a macro- 
scopic quantum self-trapping behavior of polaritonic mo- 
des takes place, first, as a result of nonlinear interaction 
of macroscopically large number of polaritons from dif- 
ferent branches. Second, the MQST regime occurs as a 
result of the "external force" action (see (fT8|h (|23| ) de- 
pending on effective difference AE c s 7^ in (fl~3|) . 

In first case, the dependence on normalized time r for 
fractional population imbalance <r is presented in Fig. 2 to 
visualize the various MQST regimes. In fact, the MQST 
effect at the phase 6 = occurs for several cases: in two 
regions 3,5 in Fig. la, in regions 3,4 in Fig. lb, and in re- 
gion 3 in Fig.lc. The curve 1 in Fig. 2 demonstrates a self- 
trapping behavior of polariton interaction for last case. 
The regimes discussed here appear for localized (periodi- 
cally) phase. A similar dynamic MQST behavior occurs 
for nonlinear interaction of two atomic BEC modes (fun- 
damental and first excited) under the coupled two-mode 
approximation - cf . [10] . 

The second kind of the MQST regime associated with 
so called it -phase MQST mode is presented in Fig.l for 
relative phase values 9 = ±7r. The gray-color regions 
marked as 2, 4 in Fig. la and as 2 in Fig.lc are respon- 
sible for that and characterize the 7r-phase pendulum 
rotations with closed loop trajectories around the value 
(?) ^ - cf. [5]. In Fig. 2 the curves 2,3 demonstrate 
the MQST dynamics discussed in respect of the regions 
taken from the Fig. la and Fig.lc. The curve 3 describes 
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the large amplitude MQST oscillations appearing here 
for the time-averaged population imbalance (?) 7^ dis- 
placement as a result of presence of zero-point energy 
difference AE cS ^ 0. 

We are emphasis here the two types of the 7r-phase 
modes occurring in the MQST regions in Fig.la,c for the 
phase values 9 = ±7r. Two various dynamical trapping 
regimes (above and beyond the stationary ^-symmetry 
breaking value of population imbalance ? s ) have been 
firstly discussed for macroscopic atomic trapping in the 
BJJ problem [5]. The modified conditions for that can 
be represented as follows (cf. [6]): 

(?) < ?(0), (31) 

for 1 type of MQST, and 

(?) > ?(0) (32) 

for 11 type of MQST. At the stationary (threshold) point 
? s we have (?) = ?(0) = ? s , i.e. the average population 
imbalance coincides with it initial value ?(0). 

In contrast with the BJJ problem in our case for 
atomic condensate 7r-phase modes one can obtain a 
new type of zero-phase MQST modes when the system 
demonstrates the self-trapping effect according to the 
conditions (|3Tj) . (f32|) . In fact, we speak about the MQST 
states appearing in the regions 3,5 (Fig. la), in regions 
3,4 (Fig. lb) and in region 3 (Fig.lc). For stationary 
points Si. 4 in Fig. lb we have (?) = ? s . 

Finally, third kind of the MQST dynamical regime 
occurs for various running phase modes with unbounded 
phase 9. The case is presented by white-color regions 1 
in Fig.l. The curve 4 in Fig. 2 characterizes the MQST 
appearing in the region 1 in Fig. lb. 

4-4 Switching dynamics for polaritons 

The MQST regimes under consideration can be achieved 
by different ways depending on the crucial parameters 
j3, A and also on the values of initial conditions for rela- 
tive phase 8(0) and fractional population imbalance ?(0). 
Practically we can vary one of them keeping others un- 
changed. For example, let us consider transition from the 
regime 5 of the Rabi oscillations in Fig. lb upward, i.e. 
to the MQST regime 1 with unbounded running phase. 
We can establish the condition for such a transition in 
this particular case in the form: 

H (? (0) , 8 (0)) > H scp , (33) 

where H (?(0), 9(0)) is the initial (conserved) energy of 
the system defined in Q12p . the energy H scp = — 4 (1 + A) 
characterizes the state that corresponds to two separa- 
trixes going through the points P\_i in Fig. lb. For in- 
stance, in the case when the parameters A, (5 and initial 
phase 9(0) are unchanged the bound energy H scp = 0.5 



5 6- 
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(b) 



Fig. 3 Time-averaged population imbalance (?) versus (a) - 
? (0) /? c and (b) - /3-parameter. The parameters for (a) are: 
? c ~ 0.2588, (5 = 1, and for (b) - ? (0) = 0.2. In all cases 
parameter A — —2 and initial phase 6 (0) = 0. 

determines a critical value of fractional relative popula- 
tion imbalance ? c w 0.2588. The MQST dynamics occurs 
in the case when initial value ?(0) > ? c . 

In Fig. 3a we plot the dependence for time-averaged 
population imbalance (?) as a function on ratio ?(0) /? c . 
The value (?) = in Fig. 3a characterizes a polariton 
interaction for the region 5 in Fig. lb. Sharp changing of 
population imbalance at ? (0) = ? c implies the "phase 
transition" to the MQST state with (?) ^ (region 1). 
By similar manner one can obtain a switching effect for 
dynamical regimes by varying of the yl-parameter. 

In Fig. 3b we establish a switching phenomenon for 
the time-averaged population imbalance (?) depending 
on /3-parameter (i.e. on kinetic energy of polaritons in 
transverse plane). The points A,B,C mark the values of 
/3-parameter for the phase portraits presented in Fig.l. 
Initial value of (?) at (3 = characterizes the aver- 
age fractional population imbalance for polaritons along 
the cavity length when the supercurrent I c = - see 
Eq. (|16p . The curve AB determines the transition from 
the MQST regions in Fig. la to the region 5 of the Joseph- 
son oscillations in Fig. lb. Then, with increasing the ab- 
solute value of /3-parameter, the regions 5 and 4 are uni- 



Josephson dynamics for coupled polariton modes under the atom-field interaction in the cavity 



9 



fied in phase plane and form a new drop- like region where 
the system is self-trapped (is not shown in Fig.l). Such 
an evolution of the phase portrait in phase-plane is ac- 
companied by sharp transition presented in Fig. 3b near 
the value of (3 ~ 1.2. Therefore, the drop-like region is 
collapsed to stationary point 5*2 in Fig.lc with further 
expanding of the /3-paramcter. The curve BC in Fig. 3b 
describes a switching effect from the MQST state with 
locked phase to the running phase MQST states deter- 
mined by the region 1 in Fig.lc. 

The switching dynamics of polaritons depending on 
the /3-parameter has an important feature in the frame- 
work of polariton condensation problem for atomic sys- 
tem. In particular, for initially non-equilibrium state of 
polaritons in the cavity with large momentum fcy , the 
variation of yl-parameter can be used in experiment to 
achieve a quasi-condensation with fen = for low-branch 
polaritons. In turn, the last parameter can be easily var- 
ied in the experiment by changing the total number iV ex 
of polaritons in the cavity (we can change, for example, 
the intensity of laser field interacting with the atoms). 
Thus, observed dynamical regimes for both fractional 
population imbalance ? and relative phase 8 allow, first, 
to estimate the fraction of condensed polaritons in the 
cavity taking into account the relevant expression (|16|) 
for the current I c and, second, to manipulate with those 
coherence properties. 

5 Conclusion 

In the paper we developed a quantum theory of Joseph- 
son dynamics for two macroscopic polariton modes (be- 
longing to upper and lower energy branches being cou- 
pled together due to Kerr-like nonlinearity of atomic 
medium) . The phenomenon can be used to examine the 
macroscopic superfluid properties of condensed polari- 
ton state as well. The analogy with nonrigid pendulum 
behavior has been carried out to understand the differ- 
ent regimes of polariton interaction. We found out the 
dispersion relations and the characteristic frequencies for 
various regimes of collective excitations in the cavity in- 
cluding the several effects, i.e. the Josephson and Rabi 
oscillations, different kinds of the MQST states for frac- 
tional relative population imbalance Some of them 
can be represented as the new one, and have not been 
considered yet. 

The realization of macroscopic tunneling and inter- 
branch interaction for polariton modes is important to 
study in experiment the relative phase coherence in the 
presence of polariton condensation (and/or quasi-con- 
densation) . 
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